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INTRODUCTION 


"Time  series  analysis  is  now  widely  used  in  asny  branches 
of  engineering,  Che  physical  sciences  and  economics .  One  important  aspect 
of  time  series  analysis  is  spectral  analysis,  which  is  concerned  with  the 
splitting  up  of  time  series  into  different  frequency  components. 
Applications  of  spectral  analysis  cower  a  wide  range  of  problems,  for 
example,  the  effect  of  wave  oscillations  on  the  vibration  of  ships  and 
the  influence  of  disturbances  or  noise  on  the  performance  of  electrical 
guidance  systems  and  chemical  reactors'*  [1]. 

This  Memorandum  provides  a  documentation  of  a  computer 
program  which  uses  Advanced  Continuous  Simulation  Language  (ACSL*)  for 
stochastic  time  series  analysis  on  the  ARL  SEC  System-10  computer.  The 
/formulae  and  terms  used  in  the  program  are  defined  in  section  2,  while 
/  section  3  provides  a  description  of  available  lag  windows  and  discusses 
/  window  closing  pnd  window  carpentry  techniques.  An  outline  of  the  computer 
program  SPECAN  (SPECtral  ANalysispls  provided  in  section  4,  and  an  example 
of  a  spectral  analysis  of  a  discrete  time  series  using  SPECAN  is  presented 
in  section  3. 

2.  PKWNIT10N  OF  POEMPLAE  AMD  TEEMS  USED  IN  SPECAN 

All  formulae  and  terms  used  in  the  program  have  been 
adapted  from  Jenkins  and  Natts  [1]. 

2 . 1  Nyqulst  frequency  <f„> 

The  Nyqulst  frequency  is  the  highest  frequency  which  can  be 
detected  with  data  sampled  st  intervals  of  &,  and  is  determined  by 


f„  -  1/24  (1) 


2.2  Lenxth  of  Pats  Record  (Tg) 

The  length  of  the  data  record  is  simply  defined  by 


Tr  -  HA  (2) 

where  H  is  the  total  ntaber  of  data  points  in  the  record. 

2.3  Test  for  White  Noise 

Situations  often  occur  in  practice  where  it  is  necessary 
to  test  whether  an  observed  time  series  could  be  regarded  as  a  product  of 
a  white  anise  process.  If  it  is  desirable  to  detect  whether  neighbouring 


e  Mitchell  A  Gauthier  Associates,  lac 


points  of  the  time  series  are  correlated,  then  the  application  of  952 
white  noise  confidence  Halts  to  the  autocorrelation  function  estimate 
acts  as  a  useful  test  for  white  noise. 


952  WHITE  NOISE  CONFIDENCE  LIMITS  s  1  2//N 


(3) 


2.4  Sample  Mean  (x) 


by 


The  sample  mean  of  the  complete  data  series  is  defined 


-  1  r 

*  “  N  Z 

N  1-1 


(4) 


2.5  Deviation  (x.) 

■  a 


The  deviation  from  the 


mean  of  each  data  point  is  defined 


x  “  X  -  x  ,  i  -  1 . ,N  .  (5) 

di  1 

2.6  Autocovariance  Function  Estimate  (acvf) 


If  the  observations  x. ,x_ . . x^  come  from  a  discrete  time 

series,  the  discrete  time  autocovsriance  function  estimate  is 


N~k 

C*x(k)  -  £  ^  <*i  -  *>  <*i+k  -  *> 


which  can  be  written 


cxx(k) 


_1_ 

N 


N-k 

I 

i-1 


where  k  •  0,1, . . . . ,N-1 . 


(6) 


(7) 


2.7  Autocorrelation  Function  Estimate  (acf) 

As  it  is  sometimes  necessary  to  compare  two  or  more  time 
series  with  possibly  different  scales  of  measurement,  estimates  of  the 
discrete  time  autocorrelation  function  arc  required  rather  than  estimates 
of  the  acvf.  Estimates  of  the  acf  may  be  obtained  by  dividing  the  above 
acvf  estimates  by  the  estimate  of  tbs  variance,  as  follows 


where  c  (k)  is  defined  by  equation  (7)  and  the  data  variance  corresponds 
to  the  Scvf  estiaate  for  taro  lag. 


Note;  For  a  stationary  process,  the  acf  is  a  function  of  the  lag  k  only, 
and  thus  is  essentially  defined  as  a  normalized  acvf. 

2.8  Smoothed  Spectral  Estimate  (SSE) 

Since  the  smoothed  spectral  estiaate  la  an  even  function 
of  frequency,  it  it  only  necessary  to  calculate  It  over  the  range 
0  ^  f  4  f  .  However,  to  preaerve  the  Fourier  transform  relationship 
between  the  sample  spectrum  and  the  sample  acvf,  it  is  necessary  to 
double  the  power  associated  with  each  frequency  in  that  range.  Hence 
the  formula  generally  used  to  calculate  the  SSE  is 


C_(f)  •  2A<c  (0)  +  2  I  c  (k)w(k)cos2irfkA  >  (9) 

“  1“  k-1  “  j 

where  0  <  f  4  f  ,  L  is  the  truncation  lag  number,  and  w(k)  is  the  weighting 
function  of  theNlag  window,  while  equation  (7)  need  only  be  computed  for 
k  t,  0. 


For  purpomes  of  computation  it  is  more  efficient  and 
convenient  to  assume  A-l,  in  which  case  all  aets  of  data  nay  be  processed 
alike.  Equation  (9)  nay  then  take  the  form 


£n(f) 


r  l-i 

2  i  c  (0)  +  2  r  c  (k)w(k)cos2*fk 

L"  k+1  “ 


} 


(10) 


where  0  i  f  <  1/2. 


Hote:  If  Adi,  then  the  estimate  obtained  by  equation  (9)  can  be  recovered 
from  equation  (10)  by  multiplying  by  A  and  plotting  the  estimate  against 
fA  instead  of  f. 

2.9  Smoothed  Spectral  Density  Estimate  (SSDE) 

If  acf’s  are  used  instead  of  acvf's,  the  smoothed  spectral 
density  estimate  la  computed  by  the  following  expression. 


)  L'1  t 

■  (f)  -2(1  +  2  £  r  (k)w(k)  cos2«fk  f 

“  l  k-1  “  j 


(11) 


where  0  4  f  4  1/2,  and  r _ (k)  is  defined  by  equation  (8). 

m 

Jcnkins_and  Watts  rai ns— nil  that,  in  order  to  obtain  a 
more  detailed  plot  of  i  (f) ,  the  f  values  should  be  calculated  at  a 
spacing  1/(2F),  where  F*?s  2  or  3  times  L.  The  final  formula  for  digital 
computation  of  the  SSDE  is  then 
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r  l-i 

<1+2  r 

C.  k-i 


r  (k)w(k)  co* 


»/k*7 

F  J 


where  t  ■  0,1 . F. 

3.  LAG  WINDOWS 


(12) 


The  lag  windows  available  in  SPECAN  are  selected  from 
those  widely  used  in  spectral  analysis  [1].  They  are  the  Tukey ,  Par ten, 
Bartlett  and  Rectangular  windows  and  are  plotted  In  Figure  1,  while  their 
properties  are  summarized  in  Table  1.  To  perform  window  closing/carpentry 
exercises  with  SPECAN,  the  logical  varlablea  LOOP  and  PTSPLT  aust  be  set 
.TRUE,  at  run  tine.  LOOP  la  defined  in  SPECAN  and  allows  access  to  the 
window  closing/carpentry  logic,  while  PTSPLT  is  an  ACSL  paraneter  [2] 
which  controls  the  ACSL  overlay  plot  facility.  Note  that  if  PTSPLT  is 
not  set  .TRUE.,  it  defaults  to  .FALSE,  such  that  only  the  results  of  the 
last  resolution  bandwidth  in  window  closing,  or  the  last  lag  window  in 
window  carpentry,  will  be  available  for  plotting.  It  can  be  seen  fron 
Figure  2  how  window  closing/carpentry  exercises  are  effected  in  the 
program . 


As  the  topic  of  window  carpentry  and  the  effects  of  window 
shape  on  smoothing  are  well  covered  by  Jenkins  and  Watts  [1),  only  an 
informative  an—iry  of  window  closing  and  window  carpentry  is  presented 
in  this  Memorandum.  However,  it  is  l^ortant  to  not*  that  the  critical 
decision  in  the  design  of  a  spectral  analysis  is  the  choice  of  window 
resolution  bandwidth  -  not  the  choice  of  window  shape  -  as  will  be  seen 
by  examining  the  exaaple  presented  in  section  S. 

3,1  Window  Closing 

The  stain  objective  of  window  closing  la  to  aid  the  physical 
insight  in  the  process  of  est lasting  and  interpreting  spectra.  It  involves 
confuting  saoothed  spectral  sat last** ,  initially  with  a  wide  resolution 
bandwidth,  and  then  using  progressively  saaller  bandwidth*,  such  that  the 
shape  of  the  spectrua  is  allowed  to  evolve.  The  initial  choice  of  a  wide 
resolution  bandwidth  will  usually  conceal  a  certain  amount  of  detail  in 
the  spectrua.  More  significant  detail  can  be  studied  as  the  bandwidth  is 
reduced  until  an  optima  condition  has  been  achieved  where  the  resolution 
bandwidth  of  the  window  is  less  than  the  smallest  significant  detail  in 
the  spectrua. 


Practical  problems  of  laterpratatlon  nay  arise  due  to 
Instability  of  the  estimates  (1],  however  it  is  possible  to  distinguish 
batsman  three  types  of  situation  which  may  occur  in  practice: 

(1)  It  is  sonatinas  possible  to  narrow  the  bandwidth  sufficiently 

to  reveal  most  of  tbs  significant  detail  without  Incurring 
instability. 


(ID 


Sonatinas  it  is  dear  that  la  no 
converging  to  a  stable  value. 


is*  is  tb*  spectrum 
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(lii)  A  situation  which  usually  falls  sonevhere  between  (i)  and  (11). 

These  situations  are  all  characterized  by  a  tendency  for 
the  estiaates  to  converge  initially,  but  then  to  diverge  due  to  instability 
before  definite  conclusions  can  be  drawn. 

When  performing  window  dosing  exercises  with  SPECAN,  the 
number  of  resolution  bandwidths,  HH,  desired  by  the  user  should  be  set 
in  the  nominated  data  file.  Aa  the  SSDE  for  each  resolution  bandwidth 
is  computed,  MS  is  reduced  by  one,  until  Mt-O,  at  which  stage  the  program 
will  terminate. 

3.2  Window  Carpentry 

Window  carpentry  relates  to  the  design  of  the  window  which 
is  to  be  used  in  the  analysis,  and  therefore  is  of  some  importance,  albeit 
secondary  to  window  closing,  as  is  shown  empirically  in  Jenkins  and  Watts 

ci]  - 

Briefly,  Jenkins  and  Watts  state  that: 

(i)  Tukey  and  Parzen  windows  give  comparable  but  superior  results 

to  the  Bartlett  window. 

(il)  Parzen  window  has  smaller  sidelobe  effects  than  the  Tukey  and 

Bartlett  windows,  but  requires  more  autocovariances  to  achieve 
a  given  bandwidth.  This  means  that  the  spectral  estimate  will 
settle  down  to  a  steady  value  more  quickly  by  using  a  Tukey 
window  rather  than  a  Parzen  window. 

(iii)  The  Rectangular  window  performs  badly  by  exhibiting  much  larger 
sidelobe  effects  than  any  of  the  Bartlett,  Tukey  or  Parzen 
windows,  and  therefore  its  nee  is  not  recoMended. 

Bote:  The  effect  of  side  lobes  is  to  permit  values  of  the  power  spectra 
at  frequencies  distant  from  f  to  make  large  contributions  to  the  bias  at 
the  frequency  f,  creating  an  effect  known  as  "leakage". 

When  performing  window  carpentry  exercises  with  SPECAN,  Mi 
must  be  set  to  unity  in  the  nominated  data  file  as  only  one  resolution 
bandwidth  is  being  examined. 

A.  SPECTRA!  ANALYSIS  PROGRAM  -  SPECAN 

The  program  is  designed  for  use  on  the  AR1  DEC  System-10 
counter  and  is  structured  so  that  the  majority  of  computations  are 
performed  in  FORTRAN  subroutines  called  by  the  main  program,  SPECAN. CSL. 

The  program  layout  therefore  reduces  the  time  required  by  the  ACSL 
translator  and  allows  independent  access  to  the  spectral  analysis  equations 
in  the  subroutines  without  the  need  to  re-translate  and  then  re-compile 
the  aain  program.  The  dialogue  for  translating  the  main  program,  compiling 
and  loading  to  obtain  the  run  version  of  SPECAN  (at  the  teletype  terminal) 
is  provided  in  Appendix  A. 


The  Fourier  Transform  Is  evaluated  In  this  program  by 
the  Jenkins  end  Wetts  finite  difference  algorithm  [1  (app.7.11,  not  the 
Fast  Fourier  Transform  algorithm. 


Mein  Pro« 


SPECAN.CSL 


The  main  progrsa,  SPECAN.CSL  (the  .CSL  extension  Is 
required  by  the  ACSL  translator).  Is  written  in  ACSL  [2]  and  therefore 
opens  the  ACSL  library  so  that  many  alternatives  for  data  output  and 
plotting  are  readily  Bade  available  without  the  need  of  any  sddltional 
dedicated  programing.  A  working  knowledge  of  ACSL  Is  assumed. 

SPECAN.CSL  acts  as  a  co-ordinating  centre  for  all  data 
and  coiqmtat  ions  by  calling  the  subroutines  GETDAT,  MULCOR  and  AUSPEC. 

The  Nyquist  frequency,  length  of  data  record,  and  the  951  white  noise 
confidence  limits  are  calculated  directly.  A  flowchart  of  SPECAN.CSL 
Is  shown  in  Figure  2. 

A. 2  Subroutine  GETDAT 

As  the  name  implies,  subroutine  GETDAT  "gets"  the  "data" 
by  reading  the  input  data  file  specified  by  the  user,  and  transfers  the 
data  deck  and  desired  lag  window  flags  into  arrays  XN'I)  and  LH(J) 
respectively.  GETDAT  also  calls  subroutine  FILTER  and  passes  XN(I)  for 
filtering  as  desired  by  the  user. 

The  else  of  all  data  arrays  is  controlled  by  a  set  of 
DIMENSION  statements  in  the  main  program.  SPECAN.CSL.  They  are  currently 
limited  in  aise  as  follows: 

Record  of  :  data  deck  -  XN(SOO) 

deviations  -  XDBV(SOO) 

autocovariances  -  CXX(SOO) 

autocorrelations  -  RZX(SOO) 

window  weights  -  w(252) 

Smoothed  Spectral  Estimates  -  C BIX (800) 

Smoothed  Spectral  Density  Estimates  -  MIX  (800) 

selected  lag  windows  -  LW(4) 

4.3  Subroutine  FILTER 

Called  by  GETDAT ,  FILTER  is  a  dummy  subroutine  which  allows 
the  user  to  incorporate  a  number  of  filter  options  as  may  be  desired.  The 
required  filter  may  than  be  selected  by  aattlng  the  parameter  P  in  the 
data  fUe  to  the  appropriate  Index.  P-1  is  nominated  as  a  default  value 
and  Indicates  that  no  filtering  of  the  data  is  required  and  that  all  the 
data  is  to  be  used. 

4.4  gubrouno-  «anjna 


MDLCOR  is  based  on  the  flowchart  MDLT1C0R  of  Jenkins  and 
*****  (app-  *3)].  Tb*  calculations  of  tho  sap  la  ittn  mi  tha  deviation* 
*  the  mss'  ta  aquations  (4)  and  (5)  respectively.  The  calculations 
«*r  st  last  or  and  tbs  acf  use  aquations  (7)  and  (»  respectively. 


...  h 


4.5  Subroutines  AUSPEC 


AUSPEC  is  based  on  the  flowchart  AUTOSPEC  of  Jenkins  and 
Watts  [1  (app.7.1)].  The  SSE  and  SSDE  calculations  use  equations  (10) 
and  (12)  respectively  and  the  selected  window  is  called  according  to 
the  parameter  LASWIN. 

4.6  Window  Subroutines 


Subroutines  SECT,  BART.  TUKEY  and  PARZEN  represent  the 
Rectangular,  Bartlett,  Tukey  and  Parzen  lag  windows  respectively.  The 
weighting  functions  for  these  windows  are  as  given  in  Table  1,  a  sumary 
of  lag  window  properties.  The  resolution  bandwidth,  degrees  of  freedom, 
and  asymptotic  variance  are  all  calculated  within  the  appropriate 
window  subroutine. 

4.7  Input  Data  Description 

The  following  parameters  are  used  as  inputs  and  are  made 
available  to  the  main  program  by  the  subroutine  GETDAT : 

(a)  N  »  total  number  of  values  in  the  data  deck. 

(b)  NCOLS  -  number  of  coltmis  of  the  data  deck. 

(c)  MN  -  the  number  of  resolution  bandvldths  to  be  used  In  windowing 

exercises.  (It  is  suggested  that  MN-3  be  used  normally  In 
window  closing  exercises  in  order  to  avoid  confusion  in  the 
overlayed  spectral  plot.  Set  Mi-1  for  window  carpentry 
exercises) . 

(d)  MAXM  •  a  constant  which  defines  the  maximum  value  of  the  truncation 

lag  number,  L.  (MAXM  is  used  in  the  switching  logic  in 
SPECAN  to  determine  whether  window  closing  or  window  carpentry 
exercises  are  to  be  performed  once  LOOP  has  been  identified 
as  .TRUE.  In  window  closing  exercises,  L  is  successively 
halved,  as  can  be  seen  in  Figure  2,  whereas  in  window 
carpentry  exercises,  L  remains  constant  and  equal  to  MAXH. 

It  is  recommended  that  the  initial  value  of  MAXM  be  chosen 
to  be  approximately  H/4,  noting  that  SPECAN  requires  MAXM 
to  be  a  multiple  of  2°  if  Ml  -  nfl,  where  n  is  an  integer). 

(e)  DELTA  -  sampling  interval. 

(f)  HP  •  the  number  of  frequency  values  at  which  the  spectral  estimate 

is  to  be  considered,  (normally  2-to-3  times  MAXM  [1]). 

(g)  P  is  a  parameter  specifying  which  filter,  if  any,  is  required. 

Default  P-1  :  mo  filtering,  all  data  is  used.  (Filters 
to  be  supplied  by  users) . 

(h)  LW( 1) ,  LW(2) ,  LW(3) ,  LW(4)  :  the  lag  window  selection  indices  In 

order  of  use.  (The  windows  are  selected  according  to  the 
parameter  LAGVLN  In  the  main  program,  as  follows: 

LACVDi  •  LW(J) ,  J-l . 4 

where  LV(J)  -  1  :  Parzen  window 
Uf(J)  -  2  :  Tukey  window 
LH(J)  -  3  :  Bartlett  window 
Ui(J)  -  4  :  Rectangular  window). 
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For  exanple:- 

(1)  10  0  0 

(ii)  4130 

(ill)  0234 


window  closing  (MN>1)  using  Parzen  window; 
no  window  carpentry, 
ie.  LAGWIN  -  LW( 1)  -  1. 

window  closing  (MN>1)  using  Rectangular 
window;  or  window  carpentry  (HN  >1)  between 
Rectangular,  Parzen  and  Bartlett  windows, 
le.  LAGWTN  -  LW( 1 )  -  4 

-  LU(2)  -  1 

-  LU(3)  -  3. 

LAGWIN  »  LW( 1)  ■  0,  is  a  run  termination 
condition  as  no  lag  window  has  been  specified 
for  first  use.  (Note  that  LAGWIN  •  LW(J)  4  0 
and  LAGWIN  “  LW(J)  >  4  are  general  termination 
conditions) . 


(i)  XN(1),  XN(2) . XN(N):  the  data  deck.  Great  care  should  be  taken 

to  avoid  data  errors  as  the  spectrum  can  be  quite  sensitive 
to  them. 

4.8  (Nitput  Data  Description 


Important  statistical  information  and  select  data  parameters 
are  automatically  printed  into  the  ACSL  print  file  PDR20.DAT,  providing  an 
echo-check  of  N,MN,MAXM, DELTA, NF,P  and  LAGWIN,  as  well  as  recording  the 
following  information: 


(a) 

XMEAK 

(b) 

CXX(l) 

(c) 

FN 

(d) 

TR 

(e) 

BWIDTH 

(f) 

DOF 

(S) 

AT 

(h) 

UWNL 

(i) 

L 

-  the  data  mean,  x,  as  calculated  by  equation  (4). 

-  the  data  variance  for  zero  lag,  c^(0) »  as  calculated  by 
equation  (6). 

•  the  Nyquist  frequency,  f^ ,  as  calculated  by  equation  (1). 

-  the  total  length  of  the  data  record,  T^,  as  calculated 
by  equation  (2). 

-  the  resolution  bandwidth,  as  defined  in  Table  1. 

■  the  degrees  of  freedom,  as  defined  in  Table  1. 

•  the  asymptotic  variance,  as  defined  in  Table  1. 

-  the  Upper  White  Noise  Limit,  corresponding  to  the  95X 
white  noise  upper  limit,  as  calculated  by  equation  (3). 

■  the  truncation  lag  number  for  the  respective  resolution 
bandwidth. 


Additional  output  may  be  selected  at  run-drive  time  by  ACSL 
commands  (2]  provide  one  or  more  combinations  of  the  following  output 
possibilities: 

(I)  data  printout  and  plot  (XT  versus  DATUM,  with  or  without 
XMEAN  over laved); 

(II)  deviation  printout  and  plot  (XTDEV  versus  DATUM); 


(ill)  acvf  estimate  and  plot  (QCXK  versus  R)  ; 

(Iv)  acf  estiaate  and  plot  (RXXK  versus  R,  with  or  without  the 

95!  white  noise  confidence  lialts  overlayed); 

(v)  lag  window  weighting  printout  and  plot  (WK  versus  K)  for 
individual  windows  or  for  window  carpentry; 

(vi)  SSE  printout  and  plot  or  log  plot  (CKXXF  versus  F) ; 

(vil)  SSDE  printout  and  plot  or  log  plot  (RBXXF  versus  F) ; 

(viii)  options  (vl)  and  (vil)  My  be  plotted  as  individual  plots  or 
overlayed  for  window  closing  using  several  bandwidths. 

The  above  options  are  the  aost  coamonly  used  [1].  However, 
the  use  of  ACSL  provides  the  user  with  a  vide  range  of  output  options  and 
nodes  of  presentation,  which  are  not  normally  available  without 
considerable  additional  dedicated  prograaaing.  Plots  which  are  requested 
are  automatically  plotted  into  the  ACSL  plot  file,  FOR22.DAT. 

A. 9  Validation  Test 


In  order  to  attain  a  reasonable  level  of  user  confidence 
In  the  progran,  a  suitable  practical  example  given  in  Jenkins  and  Watts 
[1],  namely  the  spectral  analysis  of  the  batch  data,  was  used  as  a  test 
case  to  validate  SPECAN.  The  results  obtained  from  SPECAN  agreed 
extremely  well  with  those  provided  in  Jenkins  and  Watts. 

5.  AH  EXAMPLE  OF  A  SPECTRAL  ANALYSIS  PSINC  SPECAN 


The  following  exasple  is  a  brief  spectral  analysis  of 
some  aircraft  vibration  test  data  and  is  presented  here  in  order  to 
Illustrate  the  use  end  scope  of  SPECAN.  Evaluation  and  interpretation 
of  the  results  obtained  Is  limited  to  that  which  Is  considered  necessary 
to  describe  certain  characteristics  of  SPECAN  and  ACSL. 

The  data  for  this  analysis  arose  from  acceleroweters 
responding  to  a  bonker  rocket  firing  to  port  with  200  lb  thrust  for 
50  msec.  The  acceleroweters  were  mounted  on  the  aircraft  fin  and  rudder, 
while  the  bonker  rocket  was  mounted  at  the  rudder  horn  balance.  The  data 
was  recorded  at  Intervals  of  20  msec,  and  therefore  represents  an  example 
of  a  discrete  time  series  as  values  were  only  obtained  at  specific 
intervals  of  time.  The  Input  data  file  APT. OAT  for  the  preliminary 
analysis  is  Included  In  Appendix  B. 

5.1  Preliminary  Analysis 

By  choosing  the  input  data  values  with  careful  regard  to 
section  A. 7,  it  is  possible  to  conduct  a  complete  preliminary  analysis 
with  one  run  of  SPECAN.  The  run-drive  cnMinds  for  the  preliminary 
analysis  are  listed  In  Appendix  B. 
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5.1.1  Data 

Figure  3  Illustrates  the  time  aeries  of  the  aircraft 
vibration  test  data  as  generated  by  SPECAN,  with  the  data  mean  shown 
overlayed.  As  no  obvious  trends  or  periodicities  were  revealed  and 
the  tlae  series  appears  reasonably  unifora  about  the  aean,  the  data 
is  assuaed  to  be  representative  of  a  stationary  tlae  series  and  filtering 
may  be  considered  unnecessary. 

5.1.2  Pilot  Spectrua 

An  lanediate  consequence  of  the  statlonarity  assuaptlon 
above  is  that  both  the  acvf  and  the  acf  are  functions  of  only  the  lag  [1]. 
It  therefore  makes  little  difference  which  of  the  two  functions  Is  used 

to  determine  the  spectrua,  as  the  acf  Is  simply  a  normalized  acvf,  as 

defined  In  section  2.7.  The  acf  is  used  here  as  an  Intermediate  step  In 
the  estimation  of  the  spectral  density  function,  and  also  as  a  guide  for 
designing  a  spectral  analysis. 

As  N-250  data  points,  the  initial  value  of  MAIM  was  chosen 
as  60  by  the  relationship  recommended  In  section  4.7(d),  hence  the  acf 

estimate  was  computed  up  to  L-60  lags  in  SPECAN  to  achieve  the  pilot 

spectrum  shown  in  Figure  4.  Note  that  the  spectral  density  has  been 
plotted  on  a  logarithmic  scale.  The  Tukey  lag  window  was  used,  and  NF 
was  chosen  as  ISO  In  accordance  with  section  4.7(f). 

The  bandwidth  and  80S  confidence  interval  shown  on  the 
pilot  spectrua  constitute  a  pair  of  "tolerance  limits".  The  conversion 
from  degrees  of  freedom  to  confidence  interval  is  explained  by  Jenkins 
and  Hatts  (1],  and  the  conversion  graph  from  that  reference  is  presented 
as  Figure  5. 


Important  points  to  aote  from  Figure  4  and  the  computer 
print -out  for  the  preliminary  analysis,  which  is  included  in  Appendix  V, 
are  that: 

(1)  the  AT  la  quite  large  (0.180),  which  suggests  that  some  of 

the  smaller  peaks  aay  be  spurious; 

(il)  The  DOF  is  small  (11),  such  that  the  80S  confidence  interval 

is  quite  wide ,  indicating  that  the  large  peaks  are  real,  while 
the  validity  of  the  others  is  questionable; 

(ill)  MOTH  is  quite  small  (0.0222  cps)  with  respect  to  the  Important 
detail  in  the  spectrum,  indicating  poor  stability  or  instability. 

5.1.3  Saile  Autocorrelation  Function 

The  sample  acf  for  the  aircraft  vibration  test  data  is 
plotted  In  Figure  6,  while  a  computer  print-out  of  the  preliminary  analysis 
in  Appendix  b  lists  the  first  60  lags.  The  95S  white  noise  confidence 
limits  are  shown  superimposed  over  the  sample  acf  in  Figure  7  as  a  test 
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for  white  noise.  For  the  window  closing  exercise,  the  truncation  lag 
numbers.  L,  of  8,  16  and  32  have  been  chosen,  where  MAIM  »  32;  noting 
the  discussion  in  section  4.7(d). 

5.2  Smoothed  Spectral  Density  Estimates 

The  SSDE  were  computed  while  the  Tukey  window  was  dosed 
through  the  selected  values  of  L.  They  are  plotted  together  in  Figure  8 
by  using  the  ACSL  FTSPLT  facility.  The  window  resolution  bandwidth  and 
80Z  confidence  intervals  are  shown  for  each  truncation  lag  number.  A 
computer  print-out  of  the  results  from  SPECAN  is  Included  in  Appendix  U. 
The  number  of  frequency  values.  NF,  was  chosen  as  80  in  accordance  with 
section  4.7(f).  so  that  for  this  run  of  SPECAN  the  following  input 
parameters  were  adjusted: 


(i) 

in  AVT.DAT 

:  MAIM  - 

32  .  NF  -  80 

(ID 

at  run-time 

:  FTSPLT 

-  .T. 

LOOP 

-  .T. 

CINT 

-  0.02. 

Note:  CINT  is  the  ACSL  communication  interval  parameter  [2]  and  is  set 
to  the  default  value  of  1.0.  SPECAN  is  written  such  that  if  the  sampling 
interval  of  the  data,  DELTA,  is  not  equal  to  unity,  than  CINT  must  be  set 
at  run-time  to  be  equal  to  DELTA,  otherwise  the  run  will  be  aborted  and 
an  error  message  will  be  displayed,  as  shown  in  Figure  2. 

5.3  Interpretation  of  the  Smoothed  Spectral  Penalty  Estimates 

5.3.1  Fidelity 

From  an  inspection  of  Figure  8,  it  is  clear  that  the 
resolution  bandwidth  corresponding  to  the  estimate  for  L“8  is  too  wide 
to  reveal  all  the  detail  in  the  spectra,  while  the  estimate  for  L«16 
tends  to  under-estimate  the  peaks.  If  the  width  of  a  peak  is  measured 
by  the  distance  between  the  half  power  points  [1],  it  can  be  seen  in 
Figure  8  that  the  resolution  bandwidth  for  MAIM  -  L  a  32  is  slightly  less 
than  the  width  of  the  narrowest  peak,  and  therefore  high  fidelity  is 
implied. 

5.3.2  Stability 

The  variance  of  the  estimates  is  an  indication  of  stability. 
The  variance  decreases  with  increasing  DOF  (Figure  5)  while  the  DOF  in 
turn  Increases  with  the  mater  of  samples,  N  (Table  1).  On  the  other  hand, 
DOF  decreases  with  increasing  lag  number,  L.  This  contrasts  with  the 
resolution  bandwidth,  in  section  5.3.1  above,  which  decreases  with 
increasing  L.  Thus  the  choice  of  the  truncation  lag  number  reflects  a 
compromise  between  Fidelity  and  Stability. 

Note:  As  the  SSDE  for  L>32  is  relatively  smooth  and  exhibits  high  fidelity 
aad  good  stability,  there  appears  to  be  no  benefit  la  closing  the  Tukey 
window  any  further. 
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5.3.3  Resolution 


An  indication  of  good  resolution  is  the  fact  that  the 
resolution  bandwidth  for  L>32,  using  the  Tukey  window,  is  slightly  less 
than  the  width  of  the  narrowest  Important  peak  in  the  spectrua.  In 
general,  high  fidelity  implies  good  resolution. 

5.3.4  Bias 

Bias  is  a  measure  of  the  offset  distortion  of  the  smoothed 
spectral  estimate  from  the  true  spectrua.  The  total  error  in  the  spectral 
estimate  is  a  combination  of  the  variance  of  the  estimate  and  the  bias 
error.  While  variance  may  be  decreased  (and  stability  increased)  by 
increasing  the  resolution  bandwidth  (decreasing  L) ,  this  has  the  effect 
of  increasing  bias.  A  compromise  is  therefore  required,  keeping  in  mind 
also  the  requirement  for  good  resolution  (and  fidelity)  discussed  above. 

The  topic  of  bias  is  aore  fully  covered  in  Jenkins  and  Watts  [1]. 

5.3.5  Aliasing 

As  the  topic  of  aliasing  with  respect  to  spectral  analysis 
is  adequately  discussed  by  Jenkins  and  Watts  [1]  and  other  authors,  such 
as  Kanasewich  [3],  and  Bendat  and  Pieraol  [4],  it  will  be  suffice  to  say 
that  inspection  of  Figure  8  does  not  reveal  any  detail  to  suggest  the 
presence  or  possibility  of  aliasing. 

5.3.6  Frequency  Limit 

The  highest  frequency  which  nay  be  reliably  detected 
corresponds  to  the  Hyquist  Frequency. 

5.4  Example  of  Window  Carpentry 

Recall  from  section  3.2  that  the  Tukey  and  Parsen  windows 
are  comparable.  It  can  be  shown  in  a  similar  way  to  the  following  comparison, 
that  these  windows  have  very  similar  performances.  For  this  reason,  and 
because  results  from  the  Tukey  window  have  already  been  illustrated  in 
Figure  8,  the  Tukey  window  has  been  omitted  from  the  following  window 
carpentry  example. 

Using  the  data  of  the  aircraft  vibration  test,  the  SSDE's 
were  computed  for  the  Rectangular,  Bartlett  and  Parsen  windows,  and  are 
plotted  together  in  Figure  9.  It  is  readily  seen  that  the  Rectangular 
window  performs  badly  by  exhibiting  mich  greater  Instability  characteristics 
than  either  the  Bartlett  or  Parsen  windows,  which  display  reasonable  shapes 
and  have  similar  performances. 

When  comparing  the  results  obtained  in  window  closing 
(Figure  8)  against  those  in  window  carpentry  (Figure  9) ,  it  is  necessary 
to  a^hasize  that  it  is  the  choice  of  window  resolution  bandwidth  -  not 
the  choice  of  window  shape  -  that  is  the  i^ortant  question  in  the  design 
of  a  spectral  analysis.  A  co^ratsr  printout  containing  the  i^ortant 
information  associated  with  Figure  9  is  Included  in  Appendix  B. 


A  spectral  analysis  of  data  provided  from  an  aircraft 
vibration  test  has  been  conducted  as  an  example  in  the  use  of  SPECAN. 

Filtering  was  considered  unnecessary  as  no  trends  or  periodicities  were 
detected  in  the  dsta.  The  pilot  spectrum,  obtained  as  a  preliminary 
analysis  using  a  truncation  lag  number  of  60,  exhibited  high  fidelity 
but  displayed  signs  of  instability.  The  sample  acf  was  examined  and 
new  truncation  points  selected  such  that  MAXM  ■  32  was  used.  The  Tukey 
lag  window  was  then  closed  over  the  SSDE,  resulting  in  the  estimate 
corresponding  to  L-32  being  identified  as  possessing  high  fidelity,  good 
resolution,  good  stability,  low  variance  and  reasonably  small  bias. 

Window  carpentry  was  illustrated  and  the  importance  of 
choosing  the  correct  resolution  bandwidth  (truncation  lag  number),  rather 
than  the  window  shape,  was  emphasized .  The  Rectangular  window  was  shown 
to  be  unsuitable  for  this  anal y a is  because  of  its  instability  characteristics. 

6.  CONCLUSION 

A  program  in  ACSL,  which  has  been  written  to  perform  the 
analysis  of  stochastic  time  aeries  on  the  ARL  DEC  System-10  computer,  has 
been  documented  in  this  Memorandum.  The  program  has  been  based  on  the 
text  of  Jenkins  and  Watts,  and  has  been  previously  validated  by  using  an 
example  provided  therein.  An  example  of  a  spectral  analysis  using  the 
program  has  been  presented.  The  use  of  ACSL  requires  little  programing 
effort  and  provides  the  user  with  flexibility  in  selecting  various  output 
options  with  different  modes  of  presentation. 
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APPENDIX  A 


DIALOGUE  FOR  LOADING  SPECAN  AT  TELETYPE  TERMINAL 


To  obtain  a  run  version  of  SPECAN,  the  procedure  illustrated 
below  should  be  followed  when  using  the  teletype  terminal .  Note  that  the 
uaer  subroutines  are  grouped  together  in  the  file  SPECTL. SUb. 


. RU  nuS:ACSL 


•SPECAN 

ACSL  NOBEL  PILE:  SPECAN.CSL 
LIST  ON  PILE:  SPECAN. LSI 
CONPILE  FRON:  SPECAN. FOR 
LIST  OPTNS  USED:  ET 

END  OF  EXECUTION 

CPU  TINE:  16.20  ELAPSED  TINE:  56.82 
EXIT 

.CON  SPECAN 
FORTRAN:  SPECAN 
NAIN. 

-Rloc« 

IZSlNL 

Z1BERV 

.CON  SPECTL.SUB 

FORTRAN:  SPECTL 

GETBAT 

FILTER 

NULCOP 

AUSPEC 

RECT 

BART 

TUNE* 

PAR ZEN 

.LOAD  SPECAN, SPECTL. NUS:CAPPlT,NUS:ACSLII.'SEARCH 
LIN*:  Loading 

exit 

.SAVE 

SPECA*  saved 


I 

< _ 

1 

j. 

APPEND LX  B 

'  «  SPECTRAL  ANALYSIS  OF  AIRCRAFT  VIBRATION  TEST  DATA: 


1% 

INPUT  AND  OUTPUT  BATA 

r- 

SECTION 

B.l 

Data  file  AVT.DAT  for  preliminary  analysis. 

"V 

,T 

B.2 

Run-drive  cosmands  for  preliminary  analysis. 

v"- 

B.3 

Results  of  preliminary  analysis. 

t 

v; 

B.4 

Results  of  window  closing. 

B.5 

Results  of  window  carpentry. 

B.l 


B.  1  Data  file  AVT.DAT  for  preliAdnary 


250 

X 

NUNBER  OF 

BATA 

POINTS 

to 

: 

NUAIER  OF 

COLUANS  IN  BATA  TABLE 

1 

t 

NUNIER  OF 

BANBVIDTHS 

60 

MAXimifl  VALUE  OF  TRUNCATION  LA6  NUNBER 

0.02 

SANPLIN6  ; 

INTERVAL 

150 

NUABER  OF 

FREOUENCT 

POINTS 

t 

FILTER  TIRE 

2  0  ( 

)  0  s 

LAS  UINDOVS  DESIRED 

(1ST 

2ND..  3RD.,  ‘ 

DATA 

DECK: 

29.0 

30.0 

32.0 

21.0 

2.3 

6.5 

25.0 

38.0 

46.0 

43.0 

26.5 

28.5 

32.0 

43.0 

40.0 

3’. 5 

33.0 

27.0 

20.0 

14.5 

9.5 

16.5 

19.5 

28.0 

27.5 

27.0 

27.5 

32.3 

35.3 

39.3 

39.0 

35.0 

30.5 

31.0 

30.5 

30.0 

26.0 

26.5 

28.3 

29.5 

33.5 

31.0 

28.0 

28.0 

28.5 

29.0 

20.5 

28.0 

26.5 

27.0 

28.5 

30.0 

30.5 

30.5 

29.5 

28.0 

32.0 

35.0 

36.5 

33.0 

27.0 

24.0 

23.0 

26.0 

28.3 

23.3 

22.3 

23.3 

24.0 

28.5 

33.5 

29.0 

29.5 

30.0 

31 .0 

28.0 

26.0 

28.5 

33.3 

30.0 

26.0 

23.0 

24.3 

26.3 

25.5 

23.0 

30.0 

33.0 

30.3 

28.0 

26.5 

25.0 

28.0 

35.0 

28.0 

23.0 

18.0 

23.0 

29.3 

29.0 

29.0 

28.3 

29.5 

31.0 

26.0 

26.3 

28.0 

28.0 

28.0 

27.0 

25.5 

24.5 

25.3 

27.0 

30.0 

29.3 

29.0 

29.0 

28.5 

30.5 

29.0 

26.3 

25.5 

24.0 

29.0 

29.0 

28.3 

26.3 

25.5 

25.5 

25.0 

25.0 

26.3 

30.0 

30.0 

29.0 

24.5 

26.0 

28.0 

30.0 

29.0 

27.0 

26.5 

26.0 

31.3 

29.0 

26.0 

24.0 

23.0 

27.5 

31.0 

34.0 

30.0 

27.0 

22.0 

23.0 

24.0 

24.3 

23.3 

29.0 

28.5 

27.0 

23.0 

23.3 

29.0 

32.0 

36.5 

34.0 

30.3 

27.0 

22.5 

24.0 

24.3 

25.5 

28.0 

27.5 

27.0 

27.5 

27.5 

27.5 

25.0 

23.0 

26.3 

30.0 

35.0 

26.0 

23.5 

22.0 

24.0 

30.0 

30.5 

31.0 

31.3 

29.0 

26.0 

26.5 

27.5 

29.0 

27.0 

26.0 

22.5 

29.0 

31.3 

30.0 

27.5 

25.0 

25.5 

26.3 

28.0 

30.5 

34.0 

30.0 

26.3 

27.0 

29.0 

31.0 

29.3 

28.5 

29.0 

27.0 

25.0 

26.0 

28.0 

30.0 

30.5 

29.0 

28.0 

26.0 

30.0 

30.5 

30.5 

31.3 

30.0 

26.3 

28.0 

29.5 

30.0 

31.0 

27.3 

24.0 

29.5 

29.0 

28.5 

27.0 

26.0 

27.5 

30.0 

35.0 

31.3 

28.0 

I 
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B.2  Run-drive  cot»»nds  for  prelialnrv  analyele 


RU  SPECAN 

ACSL>SET  TITLE**AVT:  PRELIMINARY  ANALYSIS' 

AC SL > PREP AR  DATUM , X  T , XME AN , K , RXXK , UUNL , LUNL , F ,R»XXF 
ACSL>SET  C1NT*0.02 
ACSL  >START 

ZFRSOPn  File  uas  not  found 

Unit1)  DSKiFOROt . DAT/ ACCESS*SEQIN/MODE*ASC 1 1 

Enter  new  file  specs.  End  tilth  an  KALT ) 

•AUT.DATt 

ACSL>PR1NT  RXXK .DATUM. XT ,F .RBXXF 
ACSL7SET  PRMPLT*.F . 

ACSL>SET  C ALPLT * . T . 

ACSL>SET  TTLCPL*.F. 

ACSIXSET  DEFPLT*.T. 

ACSL>PLOT  •XAXIS"*DATUH,XMEAN,  H1'*50.0,'LO'*0.0 
ACSOSET  BEFPLT*.F. 

ACSL  PLOT  XTAG  '»*  -  TINE  (SECS)", XT 

ACSOPLOT  *XAX1S’«F.  XTA6  *"  <CPS)  *,RBXXF,'L06',  LO *0.07,  HI  *12.0 
ACSLXSET  X1NCPL*».0 
ACSL7SET  T1NCPL*7.0 
ACSOSET  NPCCPL*! 

ACSOPLOT  *XAXlS-*r  .  ITA6'«*  aA6>  *,  XHI'*A0.0,RXXK,'10'*-«.A,'HI'*1 .0 
ACSIXSET  BCFPLT*.T. 


i 

i 

t 

i 

i 


i 


ACSL>PL0T  'XHI  *«C.C.UUNl.  L0  *-0.A,'HJ'«1.0 
ACSDPL0T  XHI  >A0.0.lIML.  'L0  ■-0.4,'HI'*T  .0 
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ACSLXSET  ®EFPIT*.F. 

ACSLXSET  UNCPl*.F. 

ACSLXSET  STflCPL*.!. 

ACSLXPLOT  'XH!'*40.0,  XTaB  *“  (LAB)  * ,RXXK,  LO  *-0.4, 'HI "*l .0, 'CHAR'*“<“ 

ACSL>ST0P 

STOP 

END  OF  EXECUTION 

CPU  TINE:  24.73  ELAPSED  TINE:  11:21.38 
EXIT 


trnaf. 
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B.3  Results  of  preliminary  analysis 

Ac  abbreviated  version  of  the  results  froa  the  preliminary 
analysis  of  the  aircraft  vibration  test  data  is  presented  below.  The 
relevant  ACSL  output  file  is  FOR20.DAT. 


SET  T1TL£*-AVT:  PRELIMINARY  ANALYSIS" 

PREPAR  DATU«,XT,ME*N,r.mM«<..LMll,M»xXF 

SET  CINT*0.02 

START 


NUMBER  OF  DATA  POINTS  <N):  2S0 
SAMPLING  INTERVAL  (DELTA) :  .02 

MAXIMUM  VALUE  OP  THE  TRUNCATION  LAG  NUnDER  tNAXN):  60 
NUMBER  OP  FREQUENCY  POINTS  <NF>:  150 
FILTER  TYPE  SELECTED  (P)s  1 


NY0U1ST  FREQUENCY  (FN)  FOR  TRUE 
RECORD  LENGTH  (TR)  OF  DATA:  3 
DATA  REAM  (XMEANU  28.1 

DELTA:  25 
.000 

.0  CPS 

DATA 

VARIANCE  IE.  ZERO  LAG  (CXX(O)):  21. 

451544 

431  UHITE  NOISE  CONFIDENCE  EMITS  ON  ACF» 

NUMBER  OF  IANDQIOTHS  <  BIN »  :  1 

LAG  HINDOO  SELECTED  (LA6UIN):  2 

BANDWIDTH  ( BWIDTH ) :  .0222  CPS 

DE6REE5  OF  FREEDOM  (DOF):  11.1 

ASYMPTOTIC  VARIANCE  <AV):  .1800 

TRUNCATION  POINT  (U:  40 

PRINT  RXXK. DATUM, XT.F.RIXXF 

♦/- .1245 

LINE 

RXXK 

8ATUH 

XT 

F 

RBXXF 

0 

1.0000000 

0.0200000 

24.000000 

0. 

1.4823040 

1 

0.4774818 

0.0400000 

30.000000 

0.0033333 

1.0485440 

2 

0.1 7*1034 

0.0400000 

32.000000 

0.0044447 

1.4451210 

1 

-0.1813204 

0.0800000 

21.000000 

0.0100000 

1. 3338740 

4 

-0.2143142 

0.1000000 

2.3000000 

0.0133333 

0.4742279 

5 

-0. 08224*4 

0.1200000 

4.3000000 

0.0144447 

0.4482264 

28 

0.0141735 

0.5800000 

35.300000 

0.0433333 

0.7140836 

24 

—4.31 3E -04 

0.4000000 

34.300000 

0.0444447 

0.7140101 

30 

0.0030183 

0.4200000 

34.000000 

0.1000000 

0.7445344 

31 

0.0312440 

0.4400000 

35.000000 

0.1033333 

1 .0529400 

32 

0.0433422 

0.4400000 

30.300000 

0.1044447 

1.4207240 

40 

0.0024028 

1.2200000 

27.000000 

0.2000000 

2.4568150 

ISO 

0.0024028 

3.0240000 

31.000000 

0.3000001 

*. 2243352 

248  C. 0024028  4.4800000  31.300000  0.3000001  0.22*3552 

244  0.002*028  3.0000000  28.000000  0.3000001  0.22<i3552 


! 

< 


B.5 

i 


SET  PR*PLT=.F. 

SET  tuni*.l. 

SET  TTLCPL=.F. 

SET  DEFPL T  = .  T . 

PLOT  "XAXIS"=DATUH.XH£AH, 'HI  =50.0, 'LO  =0.0 
SET  DEFPL T  =  .F . 

PLOT  'XTAG  =  “  -  TIAE  (SECS)*, XT 

PLOT  “XAX IS"=F ,  XTAG  =  “  (CPS)  ",*BXXF,'LOG'.'L0'=0.O7,'HI  =12.0 
SET  X1HCPL*6.0 
SET  TIHCPL=7.0 

;■  SET  HPCCPL=1 

PLOT  -XAXIS"=X.  XTAG  =*  (LAG)  'XHI ' =40.0, RXXK ,  LO  =-0.4,'HI  =1 .0 
SET  DEFPLT=.T. 

"  PLOT  'XHI  =40.0. UUNl.  LO  =-0.4,  HI  =  1  .0 

PLOT  XHI  =40.0.LDHL.  LO  =-0.4,  HI  =1 .0 
SET  DEFPlT=.F. 

(  SET  LIHCPL'.F. 

SET  STHCPL=.T. 

PLOT  'XHI  *40.0,  XTA6  =■  (LAG)  * ,«XXK ,  L0'*-0.4,  HI ' =1 .0 , 'CHAR  ■  »-<• 
I  STOP 


I 


B.4  Results  of  window  closlm 


The  results  of  the  window  closing  exercise  In  obtelnlng  the 
aaootbed  spectrsl  density  estimates  Is  presented  below.  The  relevant 
ACSL  output  file  Is  FOR20.DAT.  Note  that  output  of  the  co-ordinates  F  and 
RBXXF  has  not  been  requested  by  either  an  OUTPUT  or  PRINT  command  and  hence 
they  have  been  suppressed. 


SET  T1TLE=*AVT:  WINDOW  CLOSING” 

PREPAR  f.RBXXf 
SET  CINT=0.02 
SET  PTSPLTs , T . 

SET  LOOP* . T . 

SET  IST0P=E' 

START 

NUMBER  Of  SATA  POINTS  CN):  250 
SANPLING  INTERVAL  (DELTA):  .02 

MAXIMUM  VALUE  Of  THE  TRUNCATION  LAG  NUMBER  CHAIN):  32 
NUMBER  Of  FREOUENCT  POINTS  CNF):  BO 
fILTER  TTPE  SELECTED  CP):  1 

NTBUIST  FREOUENCT  CFN)  FOR  TRUE  DELTA:  25.0  CPS 
RECORD  LENGTH  HR)  Of  SATA:  5.000 
data  mean  (xmean):  28.1 

DATA  VARIANCE  IE.  ZERO  LAG  (CXXC0) ) :  21.4515RV 

?5X  UHITE  NOISE  CONFIDENCE  LIMITS  ON  ACF :  ♦/-.12A5 

NUNBER  of  BANDUIDTHS  CNN):  3 
LAG  UIND0V  SELECTED  CLA6UIN ) :  2 
BANDWIDTH  CBWIDTH):  .04)7  CPS 
DEGREES  Of  FREEDOM  (DOF):  20.8 
ASTMPT0TIC  VARIANCE  (AVI:  .0?60 
TRUNCATION  POINT  (L):  32 

NUNBER  Of  BANDWIDTHS  (NN):  2 
LAG  WINDOW  SELECTED  CLA6WIN ) :  2 
BANDWIDTH  (BwIDTH) :  .0833  CPS 
DEGREES  0r  FREEDOM  (DOF):  41.7 
ASTMPTCTIC  VARIANCE  (AV):  .0480 
TRUNCATION  POINT  (L):  14 

NUMBER  OF  BanBwIDTHS  CMN):  l 
LAG  WINDOW  SELECTED  CLASUIN):  2 
BANDWIDTH  (BUIBTH):  .1647  CPS 
DEGREES  0'  FREEDOM  (DOF):  83.3 
ASTHPTOTK  VARIANCE  (AV):  .0240 
TRUNCATION  POINT  <i>;  g 

SET  PRNPLT=.f. 

SET  TTLCF-* . 1 . 

SET  CALPL**.”. 

PLOT  ”»A)IE”=f.  JTAG  *”  (CPS)  \RDXXF,aOG  ,'10  *0.07.  HI  *12.0 
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B.5  Results  of  window  carpentry 

The  results  of  the  window  carpentry  example  discussed  In 
section  5.4  are  presented  below.  Comparison  with  the  Tukey  window  can  be 
bandied  by  considering  the  results  in  Appendix  B.4  with  respect  to  those 
listed  here.  The  relevant  ACSL  output  file  is  F0R2#.DAT. 


set  title*"avt:  window  carpentry" 
prepar  f.rbxxf 
SET  CINT=0.02 
SET  LOOP*. T. 

SET  FTSPLT=.T. 

SET  IST0P=81 
START 

NUMBER  OF  DATA  POINTS  <N>:  250 
SAMPLING  INTERVAL  (DELTA):  .02 

MAXIMUM  VALUE  OF  THE  TRUNCATION  LAG  NUMBER  (NAXM ) :  32 

NUMBER  OF  FREQUENCY  POINTS  <NF>:  80 

FILTER  TYPE  SELECTED  (P):  1 

NYQUISr  FREQUENCY  (FN)  FOR  TRUE  DELTA:  25.0  CPS 

RECORD  LENGTH  (TR)  OF  DATA:  5.000 

DATA  MEAN  ( XMEAN ) :  28.1 

DATA  VARIANCE  IE.  ZERO  LAG  (CXX(OI):  21.4515?? 

751  UNITE  NOISE  CONFIDENCE  LIMITS  ON  ACF:  ♦/-.  1265 

NUMBER  OF  BANDVIDTHS  <HN):  T 
LA6  WINDOW  SELECTED  (LAGU1N):  I 
BANDWIDTH  (BUIDTH) :  .0581  CPS 
DE6REES  OF  FREEDOM  (DOF):  2?.0 
ASYMPTOTIC  VARIANCE  (AV):  .0470 
TRUNCATION  POINT  (L>:  32 

NUMBER  OF  BAMDWIDTHS  (MM):  1 
LA6  WINDOW  SELECTED  (LA6UIN):  3 
BANDWIDTH  (BWIDTH) :  .044?  CPS 
DE6REES  OF  FREEDOM  (DOF):  23.4 
ASYMPTOTIC  VARIANCE  (AV):  .0853 
TRUNCATION  POINT  (l):  32 

NUMBER  OF  BANDVIDTHS  (MN):  1 
LA6  WINDOW  SELECTED  (LA6UIN):  4 
BANDWIDTH  (BWIDTH):  .0154  CPS 
DEGREES  OF  FREEDOM  (DOF):  7.8 
ASYMPTOTIC  VARIANCE  (AV):  .2540 
TRUNCATION  POINT  (L):  32 

SET  PRNPLTs.F. 

SET  TTLCPL*.F. 

SET  CALPlT=.T. 

PLOT  "XAX 1 5"*F ,  XTAG  *'  (CPS)  ",RBXXF.  L06  ,'L0'«O.O7.  HI  *12.0 
STOP 


WINDOW  LACWIN  WEIGHTING  FUNCTION 


rXG.  4: 


Pilot  spactrua  for  th«  aircraft  vibration  taat  data 
uaine  tba  Tukay  lag  window. 
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FZG.  St  Coawtloa  plot  to  ofctaia  confidanca  intarvala  for 
tlM  aooothad  spoctral  dauity  aatioata  according 
to  tfea  nuobar  of  dagraaa  of  fraadoo,  for  (l-«)  ■ 
«.«0,  «.»5,  O.SS. 


Window  closing  to  obtain  snoothad  spectral  density 
astiaatas  tor  tha  aircraft  vibration  test  data, 
using  tbs  Tukey  lag  window. 


Smoothed  spectral  density  estimates  of  the  aircraft 
vibration  test  data,  illustrating  window  carpentry 
for  L*32. 
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